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Silicon nanocrystals with diameters between 1 and 3 nm and surfaces passivated by chlorine 
or a mixture of chlorine and hydrogen were modeled using density functional theory, and their 
properties compared with those of fully hydrogenated nanocrystals. It is found that fully and 
partially chlorinated nanocrystals are stable, and have higher electron affinity, higher ionization 
energy and lower optical absorption energy threshold. As the hydrogenated silicon nanocrystals, 
chlorinated silicon nanocrystals doped with phosphorus or boron require a high activation energy to 
transfer an electron or hole, respectively, to undoped silicon nanocrystals. The electronic levels of 
surface dangling bonds are similar for both types of surface passivation, although in the chlorinated 
silicon nanocrystals some fall off the narrower gap. 
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I. INTRODUCTION 

The manipulation of the silicon surface and its abil- 
ity to interact with molecules and radicals is gaining 
importance in view of the use of silicon nanostructures 
in hybrid inorganic-organic colloids and other functional 
materials. With a large surface-to-volume ratio, free- 
standing silicon nanocrystals (NCs) are ideal to explore 
the surface functionality. They can be obtained by ul- 
trasonic dispersion of porous silicon, 11^ liquid phase syn- 
thesis by reduction of 81014,12] or plasma processes. [2-01 
Nonthermal plasma synthesis is an efficient method for 
production of particles of mono-disperse sizes and lists, 
amongst other advantages, suppressed particle coagula- 
tion and selective heating of particles through energetic 
surface reactions. Additionally, it offers the possibility 
of dopant (P,B) incorporation during growth. |10l - [l3j 

Although silane is usually chosen as a precursor for 
plasma synthesis, SiCU has also been suggested as a 
cheaper and safer alternative. 0, Q Nanocrystals grown 
from a SiCl4/H2/Ar mixture are terminated with a 
mixture of chlorine and hydrogen, with variable frac- 
tions depending on the plasma composition and reactor 
pressure, ill, l8| 

The fraction of surface CI on silicon nanocrystals grown 
by this or other methods can also be increased by CI2 
plasma etching, treatment with a solution of PCI5 on 
chlorobenzene or with di-, tri-, and tetrachlorosilane 
gases, procedures already in use for industrial process- 
ing of silicon single crystal substrates. This leads to 
the formation of mono-, di- and trichloride terminations 
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at low temperatures (< 400°C).[T^ Monochloride is the 
most stable, remaining at higher temperature. The ad- 
sorption of chlorine and SiCl„ on flat silicon surfaces has 
been extensively studied by theoretical methods. [l6l - [2^ 
The adsorption energy of CI2 on a reconstructed Si(lOO) 
surface was found to be 5.4 eV,[l^ giving a Si-Cl bond 
energy of about 4 eV. On Si(lll) surfaces, the Si-Cl bond 
energy is similar, the calculated values ranging between 
3.5 and 4.2 eV.[17| In both cases, the barrier for chlorine 
diffusion is about 1 eV, and desorption takes place in the 
SiCl2 form.[l^, Further, it was found that with in- 
creasing chlorine supply the structure of the chlorinated 
Si(lll) surface suddenly changes, with a transition from 
a monochloride phase to a polychloride phase. [17| The 
adsorption of chlorine induces Cl-related CI- Si bonding 
states below the top of the valence band.fis'l 

The reactive Si-Cl surface bonds are convenient for sur- 
face functionalization with alkene and amine groups. [2^ 
[23 | Although Si-Cl bonds are stronger than Si-H bonds, 
the CI atoms, with a higher affinity for electrons, can 
more easily receive an electron from the highest occupied 
molecular orbital (HOMO) of the other reactant during 
the interaction. This additional electron is partially lo- 
calized on the shallowest p orbital of the CI radical, re- 
sulting in 01" being released. Forming Cl~ in the tran- 
sition state is energetically more favorable than breaking 
Si-H bonds, leading to lower activation energy barriers for 
grafting in Si-Cl bonds, even for partial 01 coverages. [2^ 

Additionally, the presence of chlorine changes the op- 
tical and electronic properties of the material, opening 
exciting possibilities for surface-driven electronic struc- 
ture engineering. Previous electronic structure calcula- 
tions have found that chlorine-covered nanocrystals have 
a lower gap between occupied and unoccupied electron 
energy levels and higher electron affinity than hydrogen- 
covered nanocrystals. [25! [26j Thus, it is possible that par- 



2 



tial or full surface chlorination can be used to control the 
position of the electronic levels for specific applications. 

Given the interest on Cl-terminated nanocrystals, both 
for subsequent surface conversion or for electronic struc- 
ture engineering, theoretical information on the stability, 
electronic and optical properties of silicon nanocrystals 
is of great interest. Therefore, we have carried out a de- 
tailed theoretical study to compare the properties of Cl- 
terminated silicon nanocrystals with 1-3 nm of diameter 
with H-terminated nanocrystals in the same size range. 
The first-principles methodology will be described in Sec- 
tion I. The structure and energetics of perfect nanocrys- 
tals with C1-, H- and mixed terminations will be consid- 
ered in Section II, and their electronic and optical prop- 
erties will be given in Section III. Section IV is dedicated 
to doped and defective CI- and H- terminated clusters. 
Finally, Section V discusses the relevance of the results. 



II. METHODOLOGY 

The electronic structure of the nanocrystals was ana- 
lyzed using first-principles calculations based on density 
functional theory, with a pseudopotential approach, as 
implemented in the AIMPRO code. QTi,i2S^ The local den- 
sity approximation [, 2 ij was used for the exchange and 
correlation energy. Core electrons were accounted for by 
using the pseudopotentials of Hartwigsen et a/..[29j 

Kohn-Sham orbitals were expanded on a localized ba- 
sis set consisting of atom-centered Cartesian Gaussian or- 
bitals with angular momentum up to ? = 2, as described 
in Ref. [sO]- For the core silicon atoms, we used a con- 
tracted basis set with 13 functions per atom (44G*), in- 
cluding a polarization function with 1 = 2, optimized for 
bulk silicon. A basis of the same size, optimized for SiH4, 
was used for hydrogen. For chlorine, an uncontracted ba- 
sis set with four I = and twelve I = 1 functions per atom 
was used. Convergence tests for silane, tetrachlorosilane 
and SigTHyg/SisTClye nanocrystals show that these bases 
offer an excellent compromise between accuracy and com- 
putational effort, specially for large nanocrystal diam- 
eters, where the electronic structure becomes increas- 
ingly bulk-like. For the worst case, the SiH4 and SiC14 
molecules, Si-H and Si-Cl bond lengths (TableU) are con- 
verged respectively within 0.006 and 0.015 A, bond en- 
ergies are converged within 0.1 and 0.5 eV, respectively, 
and the Kohn-Sham gaps are converged with 0.22 and 
0.09 eV, respectively. They are also in good agreement 
with previous LDA calculations. 

Total energy calculations were performed using finite 
real space boundary conditions. The optical absorption 
cross section was calculated using periodic boundary con- 
ditions, ensuring a minimum distance of at least 10 A 
between replicas of the system. In this case, the charge 
density was expanded in a plane wave basis set with an 
energy cutoff of 350 Ry. 

The equilibrium geometry of the nanocrystals was 
found by a relaxation of all the atomic coordinates using 



TABLE I: Bond length, bond enthalpy, Kohn-Sham gap 
{AEks ~ Sluks — EnoKs) and vertical excitation energy 
of the SiH4 and SiCU molecules. 



property Isi-x (A) AHb (eV) A^ks Ag* 



SiH4 


This work 


1.49 


3.8 


7.93 






Prev. calc. 


1.50° 


3.5' 


7.93 


8.76%9.26-' 




Exp. 


1.48'' 


3.3,"= 3.2'' 






SiCU 


This work 


2.03 


4.3 


6.37 


6.67 




Prev. calc. 








9.14» 




Exp. 


2.02'' 


4.0'' 







''From Ref. 

"All-electron LDA calculation from ReL [sjl 
''LDA calculation from Ref. 335 
<=From ReL[3j 

''From the heats of formation in Ref. |35|l 
"^GGA-PBE calculation from ReL [3^ 
^B3LYP calculation from ReLjH 
^Discrete variational Xa calculation |38ll 



TABLE II: Atomic composition, diameter and symmetry of 
the Si„Xm nanocrystals studied, where X g{H,C1}. 



n 


m 


d (nm) 


sym. 


35 


36 


1.1 


Td 


87 


76 


1.5 


Td 


244 


144 


2.1 


Td 


275 


172 


2.2 


Td 


286 


170 


2.2 


Dsd 


377 


196 


2.4 


Td 


513 


252 


2.7 


Td 


717 


300 


3.0 


Td 



a conjugate gradient algorithm. 

The optical absorption cross-section was calculated 
in the long-wavelength dipole approximation using the 
Kohn-Sham eigenvalues i?" and eigenvectors \tp'^). The 
matrix elements of r are evaluated using the momen- 
tum operator plus the commutator of the non-local part 
of the pseudopotential. [31] The Brillouin Zone sampling 
was restricted to the F point. The electronic temper- 
ature used as parameter in the Fermi-Dirac distribution 
was 0.1 eV/ks, where ks is the Boltzmann constant, and 
the Gaussian broadening used was 0.05 eV. 

III. STRUCTURE AND ENERGETICS 

A. Structure 

The nanocrystals used in this study were obtained by 
cutting an approximately spherical core out of a perfect 
silicon crystal and passivating the surface dangling bonds 
with CI or H atoms. The cutoff diameter can be esti- 
mated as d = [3n/(47r)]^/'^ao, where n is the number of 
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silicon atoms and oq is the calculated lattice parameter of 
bulk silicon (5.39 A). The surface silicon atoms were four- 
fold coordinated and had mono- or di-hydride/chloride 
termination. The number of Si and CI or H atoms in 
each nanocrystal is given in Table |TT1 All the nanocrys- 
tals were centered at an atomic site, with exception of the 
Si286^i70 nanocrystals, which were centered at a bond- 
center. 

After atomic relaxation, all silicon atoms remain four- 
fold coordinated, and the lengths and angles of the Si-Si 
bonds are close to those of the bulk crystal, specially 
those at largest distances from the surface [Fig. [T]-a) ,d) 
and g)]. Both CI- and H-terminated nanocrystals main- 
tain a marked crystalline character, characterized by a 
discrete radial pair distribution function relative to the 
nanocrystal center. However, there are quantitative dif- 
ferences between CI- and H-covered nanocrystals. As 
highlighted in Fig.[l] the bondlength distribution is much 
broader for the chlorinated nanocrystals. For example, 
for d = 1.5 nm, the average Si-Si bondlength of the 
CI- and H-terminated nanocrystals deviates only —0.01 
and -1-0. 02A, respectively, from the calculated bulk Si- 
Si bondlength (2.34 A), but the standard deviation is 
four times larger for the latter (Table IHI[) . Curiously, 
the bond angles are closer to the bulk value for the chlo- 
rinated nanocrystal, since the effective radius of the CI 
atom is closer to that of the silicon atom. 

The Si-H surface bonds are elongated about 0.2 Awith 
respect to their length in the SiH4 molecule (1.49 A, to 
be compared with the experimental value 1.48 A(39l|). 
However, the Si-Cl bondlengths are very close to those 
of SiCl4 (2.03 A, to be compared with the experimental 
value of 2.02 A^Sg). 

We now analyze the structure of nanocrystals with a 
mixed Cl/H surface. The fraction x of CI atoms was var- 
ied between and 1 {x eO.25, 0.50, 0.75), where x is 
the ratio between the number of CI atoms and the total 
number of CI and H atoms. For each x, the bondlengths 
and angles were averaged over 24 randomly generated 
samples (Table HIH) . We notice that for the three inter- 
mediate x fractions the distribution of Si-Si bondlengths 
is more narrow than for a; = 1 (Cl-covered nanocrystal), 
whereas the bond angle distribution is more narrow than 
for a; = (H-covered nanocrystal). 



B. Formation and reaction enthalpies 

It is important to know how the passivation with CI 
affects the stability of the nanocrystals. Although highly 
metastable structures can be prepared out of equilibrium, 
for example in non-thermal plasmas, the enthalpy of for- 
mation is an useful to characterize the stability in equilib- 
rium and estimate reaction energies. It is therefore useful 
to evaluate the enthalpy of formation of the nanocrystals 
with Cl-covered surface or with mixed CI and H surface 
passivation. 

We calculated the formation enthalpies at T = 0, de- 



TABLE III: Average Si-Si, Si-Cl and Si-H bond lengths 
{Ix~y) and angles (a), and their standard deviations, for 
a nanocrystal with diameter d — 1.5 nm, as a function 
of the surface CI coverage ratio (x). Lengths are given in 
Angstrom, and angles in degrees. The experimental bulk Si- 
Si bondlength is 2.35 A. 



X 





0.25 


0.5 


0.75 


1 


bulk Si 


'si-Si 


2.33 


2.33 


2.33 


2.34 


2.36 


2.34 


'Si-Cl 




2.10 


2.09 


2.09 


2.07 






1.72 


1.69 


1.68 


1.69 








0.0084 


0.0095 


0.014 


0.024 


0.04 




A^si-ci 




0.012 


0.044 


0.072 


0.0098 




AZsi-H 


0.39 


0.38 


0.36 


0.38 






a 


103.2 


106.0 


108.2 


109.3 


109.4 


109.5 


Aa 


20.4 


15.9 


11.2 


7.5 


4.0 





fined as the enthalpy change relative to the standard 
phases of Si (crystalline Si), CI (molecular CI2) and H 
(molecular H2) 

(1) 

where £^nc('^j ™cij ™h) is the calculated total energy of 
Si„Clm(,|Hmjj, andi?(Si), i?(Cl2) andi?(H2) are the total 
energy per atom of crystalline silicon and the total ener- 
gies of the CI2 and H2 molecules, respectively, calculated 
using the same approximations. When the formation re- 
action is exothermic, Hf is negative. 

We note that the total energies of small molecules, in 
particular H2, are not accurately calculated using the 
LDA approximation. Thus, our calculated formation en- 
thalpies of SiCU and SiH4 are underestimated: we obtain 
—7.1 and —0.2 eV, respectively, whereas the experimen- 
tal values are —6.6 and —0.4 eV[3^. However, this error 
often cancels out when calculating reaction energies. For 
example, the enthalpy change for the hydrogen replace- 
ment reaction 

SiH4 + Cl2^SiH3Cl + HCl, (2) 

which is —2.83 eV in our calculation, is only underesti- 
mated by 0.06 eV (relative to the value obtained from 
the experimental heats of formation 'SsJ). Thus, the cal- 
culated formation enthalpies can still be used to draw 
qualitative conclusions. 

Nanocrystal formation enthalpies are shown in Fig. 
a). Cl-covered clusters have lower Hf than H-covered 
clusters, typically by 1.1-1.6 eV per surface Si-H or Si-Cl 
bond. This difference is larger than the difference be- 
tween the bond energies in the SiCU and SiH4 molecules 
(Table m. It is also larger than the errors in Hf (SiCU) 
and i?/(SiH4). 

The enthalpy of formation of the clusters with mixed 
Cl/H surface follows very closely a linear interpolation 
of the endpoints x = and x = I. In analogy with the 
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FIG. 1: Calculated Si-Si bondlengths as a function of the distance of the bondcenter to the center of the nanocrystal (rec), 
and respective bond length histograms, for three nanocrystal diameters: (a-c) 1.5 nm, (d-e) 2.2 nm (Dsd) and (f-h) 3.0 nm. 
The dashed line represents the bulk Si bondlength. 



alloys, we can define a mixing enthalpy characterizing the 
deviation from linearity: 

AiJmix = Ejqcin,mci,m}i) - [£^nc("^,™ci + mH,0) (3) 

-£'nc("-, 0,TOci + m-n).] 

This is given in Fig [2]-b) for d ~ 1.5 nm. The mixing 
enthalpy AH^i^ is negative and smaller by one or two 
orders of magnitude than the enthalpy of formation. If, 
near the temperature at which the Cl- and H- atoms 
become mobile, the mixing free energy remains negative, 
this means that a binary system with Si„Clm and Si„IIm 
moieties will be unstable against the formation of a mixed 
SinClmci^mn cnscmble. 

Since the effective radius of the Cl atoms is much larger 
than that of H atoms, an additional question is whether 
steric effects prevent Cl atoms from occupying neighbor- 
ing positions, even preventing complete chlorination al- 
together. To investigate this, we calculated the enthalpy 
change associated with the hydrogen replacement reac- 
tions, 

SinClm^-.^ — iHm.jj + l + CI2 — > Si„ClmciHmH + HCl (4) 

for mci < f^, which is given by 

AHn = E^cin,mci,mn)+E{liCl)- (5) 
[ENc{n, mci - 1, ma + 1) - ^^(Cls)] 



for d — 1.5 nm nanocrystals. Two opposite situations 
were considered. Starting with SigTClHys, we first cre- 
ated an even distribution of Cl by placing each additional 
Cl atom in one of the surface sites (position r) minimizing 
the objective function 



mci — 1 

/(r)= k-rr^ (6) 

1=1 



This results in a sequence of clusters where Cl replac- 
ing for H takes place at the position further away from 
all the other Cl atoms. The enthalpy changes for this 
sequence of replacements are compared with those ob- 
tained for a concentrated Cl distribution, where each Cl 
atom is placed as close as possible to the atoms of the 
same species (thus maximizing /(r)). The results [Fig[5]- 
c)] show that there is no clear energetic preference for the 
first process, although the distribution of the enthalpies 
of replacement is smoother and narrower. Moreover, the 
enthalpies of replacement stay approximately constant 
up to 50% coverage, showing only a slight increase for 
higher x. So, there is in principle no reason why com- 
plete Cl coverage would not be attainable. 
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FIG. 3: Charge density associated with the HOKS and LUKS 
levels: (a) HOKS of Cl-covered nanocrystal (b) LUKS of Cl- 
covered nanocrystal (c) HOKS of H-covered nanocrystal (d) 
LUKS of H-covered nanocrystal. The isosurface value is 3 x 
10"^ and 8 x 10^" for HOKS and LUKS states, respectively. 
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FIG. 2: Stability of SinClmciHmH nanocrystals. (a) En- 
thalpy of formation of CI- and H- terminated nanocrystals 
as a function of the diameter, (b) Mixing enthalpy (Eq. 2)1 of 
d — 1.5 nm nanocrystals with mixed surface as a function of 
the CI fraction x = mci/(m,ci + rnn), and (c) comparison of 
the CI replacement enthalpy (Eqs. |3]- O, for an even distri- 
bution of CI (placement of each CI in the farthest position 
from the existing mci — 1 CI atoms) and for a concentrated 
distribution of CI (placement of CI in the nearest position to 
the existing mci — 1 CI atoms), for d — 1.5 nm. 



IV. ELECTRONIC AND OPTICAL 
PROPERTIES 



A. Analysis of the Kohn-Sham states 

Let us start by analyzing the electronic structure of the 
SiCl4 and SiH4 molecules, as represented by the Kohn- 
Sham eigenstates and eigenvalues. Although these quan- 
tities have only an auxiliary role in DFT, their analysis 
is useful to understand the bonding and the contribution 
of CI and H atoms to the ground state ans excited states. 
The highest occupied state of the SiCU molecule is the 
2ti state, followed closely by the 2e and 8^2 states. The 
HOKS is completely localized on the CI atoms (formed 
by CI 3p states). The 2e and 8^2 states also have a lo- 
calization of less than 10% on the Si atom. In contrast, 
the HOKS of SiH4, which is the 2t2 state, is a bonding 



state 41% localized in the Si atom. The lowest unoccu- 
pied state of the SiCU molecule is the 8ai state, followed 
1.7 eV above by the 9t2 state, whereas the LUKS of SiH4 
is 3i2- Both are partially localized on Si: 54% in the case 
of the SiCU LUKS, 65 % in the case of the SiH4 LUKS. 

As the number of Si atoms increases, the highest oc- 
cupied Kohn-Sham (HOKS) state and the LUKS start 
to resemble the bulk silicon valence and conduction band 
states, but in the case of the Cl-covered clusters the CI 3p 
character is maintained. Figure [3] depicts the charge den- 
sity associated with the HOKS (triplet) and LUKS states 
for d = 1.5 nm. The localization of those gap-edge states 
on the surface atoms is greater for the Cl-covered clus- 
ter, specially for the HOKS state (the fractional HOKS 
localizations on Cl/H are respectively 40 and 9% for the 
CI- and H-covered NCs). The contribution of the CI 3p 
atomic orbitals to the HOKS state is evident in the shape 
of the charge density isosurface near the surface of the 
nanocrystal, which resembles the SiCl4 2ti state. Simi- 
larly, near the surface the LUKS state bears some resem- 
blance to the SiCU LUKS (8ai) state. 

The calculated HOKS-LUKS gap (£:ks) of the SiH4 
and SiCl4 molecules is respectively 7.93 and 6.96 eV 
(Table |T|. Although these are not far from the experi- 
mental absorption energy thresholds, which are respec- 
tively 8.99 eVdO] and 8.84 eV,^ there are several rea- 
sons why they cannot be compared directly to experi- 
ment. Firstly, the HOKS and LUKS states of SiCU are 
2ti and 8ai, respectively, and the HOKS^-LUKS opti- 
cal transition is forbidden by symmetry; the lowest al- 
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lowed transition, 2ti 9^2, corresponds to an eigen- 
value energy difference of 8.1 eV. Moreover, the Kohn- 
Sham states change considerably in the excited state, 
and so do Coulomb, exchange and correlation interac- 
tions. Moreover, the threshold energy of the absorp- 
tion spectra of both molecules is a Rydberg transition 
(4s -J> 8i2)[43,|4li. These excitonic effects leading to Ry- 
dberg states are not described by the ground state DFT. 
To our knowledge, Rydberg transitions have not been 
resolved for undoped silicon nanocrystals. 

For silicon nanocrystals with diameters between 1 and 
3 nm, the lowest excitation energy, obtained by calcu- 
lating the difference between the total energies of each 
nanocrystal in the ground state and in the first excited 
state (at the ground state geometry), E^^ia — ~ ^G' 
differs little from £^ks (Fig. H]). This means that, upon 
excitation, the change in the electrostatic interaction en- 
ergy (resulting from the interaction between electron, 
hole, and image charges) is partially canceled by the 
change in the exchange and correlation energy. 

The excitation energy {E^) of hydrogenated nanocrys- 
tals has been previously calculated at different levels of 
theory. As illustrated in Fig. |4j our results are in good 
agreement with previous calculations with empirical po- 
tentials, and differ less than 1 eV from GW gaps. The 
reason why the LDA HOKS-LUKS bandgap is a good 
approximation for the excitation energy is clarified by 
Delerue et ai, who have proved that the differences be- 
tween the corrections to the self-energy in bulk and in 
the nanocrystal (JSi?), are nearly canceled out by the 
Coulomb interaction between electron and hole (Ec)- 
As a result, E^ — Eks + SUb, where S'Sb is the bulk 
self-energy correction, which is about 0.6-0.7 eV for the 
LDA approximation. In this work, we will assume that 
this correction is independent of the nanocrystal sur- 
face, thus justifying the comparison between CI- and EE- 
covered nanocrystals using the values directly obtained 
from first-principles. 

The minimum excitation energy is lower for the Cl- 
covered nanocrystal than for the H-covered nanocrystal. 
This follows the lower effective confinement volume for 
the HOKS and LUKS states in the Cl-covered nanocrys- 
tals. The difference is greater for the smaller diameters, 
amounting to about 1 eV for d ~ 1.5 nm. With in- 
creasing d, the gap of the Cl-covered clusters decreases 
almost linearly in this size range, whereas that of the H- 
covered nanocrystals it varies approximately with dr^'^. 
The average gap of nanocrystals with mixed CI- and H- 
coverage varies monotonically between those of the Cl- 
aud H-covered clusters of the same size (Fig. [S]). The 
variation in the gap distribution for each set of samples 
with the same d and x is not negligible, and is represented 
in Fig. [S] by the errorbars. 
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FIG. 4: Excitation energy (*) and Kohn-Sham bandgap (o) 
as a function of the nanocrystal diameter obtained from the 
present work. Data from previous calculations are shown 
for comparison: lowest excitation energy excitations cal- 
culated by solving the Bethe-Salpeter equation (GW-BSE) 
from ReflS. or obtained by time-dependent LDA (TDLDA), 
also from Refl43. calculated using empirical pseudopotentials 
without inclusion of the Coulomb interaction between elec- 
tron and hole (EP) or including this contribution (EP2), from 
Ref.43- Dotted lines are a guide to the eye. 



B. Optical spectra 

We have calculated the optical absorption cross sec- 
tion directly using the Kohn-Sham eigenstates and eigen- 
values, as described in Sec. HIl Previous theoretical work 
has shown that optical spectra can, in a good approx- 
imation, be obtained from the Kohn-Sham eigenvalues, 
since only minor charge rearrangements occur between 
ground and the low excited states, [il] As expected, the 
absorption threshold energy is lower for the Cl-covered 
cluster (Fig. [5]). However, the energies close to i?KS have 
very small or vanishing oscillator strengths. The thresh- 
old is steeper for the H-covered nanocrystal and for the 
nanocrystal with mixed surface, were the symmetry is 
broken, than for the Cl-covered nanocrystals. Overall, 
the absorption band in the 2-6 eV range is broader for 
the Cl-covered nanocrystal. 



C. Ionization energy and electron affinity 

The ionization energy was obtained from the total en- 
ergies of neutral and charged clusters, 

I{q/q + l) = E(q + l)-E{q), (7) 

where q is the charge. The results are shown in Fig. [HI 
The electron affinity of the chlorinated nanocrystals is 
higher by 2-3 eV, for the whole range of d considered, 
refiecting the higher affinity for electrons of CI. The ion- 
ization energy is also higher, but only by about 1 eV. 
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FIG. 6; Ionization energy [/(0/+)] and electron affinity 
[/( — /O)], (a) for CI- and H-covered nanocrystals, as a func- 
tion of the diameter, and (b) for d — 1.5 nm nanocrystals 
with mixed surface ( SigTClmciHrnH); ^ a- function of the CI 
fraction x = mci/{mci + mn). The absorption spectra for 
X = 0.5 was obtained from the average of ten random surface 
configurations. 



reverse is true for shallovif acceptors. However, this does 
not happen either for CI- or H-covered nanocrystals with 
d between 2 and 3 nm (Fig. [7]). This is due to the 
carrier confinement and appearance of image charges, 
which were already extensively discussed for H-covered 
nanocrystals. [45l - l49l | 



FIG. 5: Kohn-Sham bandgap (a) and calculated absorption 
spectra (b) of d = 2.2 nm silicon nanocrystals for different 
fractions of CI coverage {x = mci/{mci + mu))- The eigen- 
value gaps for fractional x were obtained by averaging over 
24 samples, and the errorbars represent the standard devia- 
tion of the results. The absorption spectra for x — 0.5 was 
obtained by averaging over 10 samples. 



V. DEFECTS 

A. Dopants 

Further, we compare the ionization energy/electron 
affinity of doped nanocrystals with that of the pristine 
(undoped) nanocrystals. This comparison is relevant 
when doping nanocrystal composites where only a small 
fraction of nanocrystals encloses one or more dopant 
atoms. In that case, the ideal is that a nanocrystal 
doped with a shallow donor (for example P) has ion- 
ization energy /(0/-I-) very close to the electron affinity 
/(— /O) of the undoped nanocrystal. Ideally, /p(0/-|-) — 
Iud{—/0), where the subscripts label the doped and un- 
doped nanocrystals should be comparable to kT, where 
T is the temperature and k Boltzmann constant. The 



B. Dangling bonds 

The position of the I{q/q + 1) levels of surface dan- 
gling bonds (DBs) relative to the gap edge states of the 
pristine nanocrystals has also been compared for both 
types of surface. We considered only dangling bonds 
on di-chloride or di- hydride surface silicon atoms i.e., 
those that in the pristine nanocrystal were attached to 
two surface terminators. In c? = 1.5 nm nanocrystals, 
monohydride DBs, although in average higher in energy 
by 0.9 eV, have similar properties. For each NC size, 
there are several nonequivalent surface Si atoms with 
di-chloride (or di-hydride) terminations where dangling 
bonds can form. The geometry and energy of each of the 
respective defects was optimized independently. 

First, we note that the (— /O) and (0/-I-) level posi- 
tions do not display a clear trend with the nanocrystal 
size. The DB levels are very dependent on the particular 
defect geometry, specially on the Cl-covered nanocrystals 
(Fig. [5]) . The main difference between CI- and H-covered 
nanocrystals is that, as the gap of the former is smaller, 
some of the DB levels fall outside the gap. That does not 
happen for the H-covered nanocrystals. 
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[/( — /O)] of doped and undoped (a) Cl-covered nanocrystals 
and (b) H-covered nanocrystals, as a function of the diame- 
ter. Shaded areas represent energies lower than the electron 
affinity or higher than the ionization energy of the undoped 
nanocrystal. The axis were inverted for easier visualization. 



VI. DISCUSSION 

Silicon NCs with chlorinated or mixed surface have 
two main potential uses: as an intermediate material for 
further surface functionalization and modification, and 
as an electronic or optoelectronic material. The energies 
and electronic properties obtained in this study may be 
useful to design or tune both types of applications. 

We have show that the formation enthalpy of the chlo- 
rinated Si-NCs is more negative than that of the hydro- 
genated Si-NCs, relative to the standard states of CI and 
H. However, this does not mean that the former will in 
general be more stable against reaction. In fact, chlo- 
rinated Si-NCs have very high electron afhnity, and can 
easily trap electrons to the LUMO state, which is par- 
tially localized on the surface CI orbitals (Fig. [3]). This 
leaves CI more susceptible to removal and substitution by 
a foreign radical. It is interesting to note that although 
mixed CI and H surfaces have a negative mixing enthalpy, 
there is no strong repulsion between nearest-neighbor CI 
atoms, specially for small coverage ratios. Thus, if the 
nanocrystals are kept neutral, there is in principle the 
possibility of engineering next-neighbor surface replace- 
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FIG. 8: Incomplete Ionization energy [/(0/+)] and electron 
affinity [/( — /O)] of surface dangling bonds on (a) Cl-covered 
nanocrystals and (b) H-covered nanocrystals, as a function of 
the diameter. Dangling bonds were formed on all symmetry- 
nonequivalent di-hydride terminations. Shaded areas repre- 
sent energies lower than the electron affinity or higher than 
the ionization energy of the undoped nanocrystal. The axis 
were inverted for easier visualization. 



ments using hydrogenated NCs with selected CI substi- 
tutions. 

Chlorine can be used to modify intentionally the elec- 
tronic and optical properties of the Si-NCs. Chlorinated 
Si-NCs also have a smaller gap between occupied and 
unoccupied electron levels. As a result, the threshold en- 
ergy for optical absorption is redshifted. The absorption 
edge can be varied by changing the CI coverage ratio. 

As both the electron affinity and ionization energy 
are greater than those of hydrogenated silicon nanocrys- 
tals, the CI surface coverage ratio can be used to tune 
the alignment between the Si-NC states and the bands 
of other materials in hetero junctions. Chlorinated Si- 
NCs doped with P [or B] have ionization energy [elec- 
tron affinity] levels quite distant from the (— /O) [(0/+)] 
levels of pristine nanocrystals with the same size. This 
also happens with hydrogenated Si-NCs, making it diffi- 
cult for a doped NC to donate free carriers to undoped 
nanocrystals. However, as the gap of the chlorinated Si- 
NCs shifted to lower energies, P(0/-|-) and B(— /O) lev- 
els are also shifted relative to the vacuum level, in com- 
parison with the hydrogenated crystals. This knowledge 
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may be useful to design heterojunctions with doped Cl- 
terminated Si-NCs as one of the components. As in the 
hydrogenated Si-NCs, dangUng bonds will act as exciton 
recombination centers. 
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